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FOREWORD 


The  study  reported  herein  was  performed  by  Dr.  Klaus  W.  Wiendieck 
of  the  Mobility  Research  Branch,  Mobility  and  Environmental  (M  and  E) 
Division,  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
as  part  of  the  vehicle  mobility  research  program  under  DA  Project  No. 

1-V-0- 14501- B52A,  "Research  in  Earth  Sciences,"  Task  1-V-0- 14501- B52A- 01, 
"Terrain  Analysis,"  under  the  sponsorship  and  guidance  of  the  Development 
Directorate,  U.  S.  Army  Materiel  Command. 

The  study  was  conducted  during  the  period  January-Becember  1966  under 
the  general  supervision  of  Messrs.  W.  J.  Turnbull,  Technical  Assistant  for 
Soils  and  Environmental  Engineering,  and  W.  G.  Shockley  and  S.  J.  Knight, 
Chief  and  Assistant  Chief,  M  and  E  Division;  and  under  the  direct  super¬ 
vision  of  Dr.  D.  R.  Freitag,  Chief,  Mobility  Research  Branch,  and  Mr.  A.  J. 
Green,  Jr.,  Chief,  Vehicle  Dynamics  Section,  M  and  E  Division.  The  report 
was  prepared  by  Dr.  Wiendieck. 

COL  John  R.  Oswalt,  Jr.,  CE  was  Director  of  WES  during  this  study, 
and  Mr.  J.  B.  Tiffany  was  Technical  Director. 
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Experimental  coefficient,  dimensionless 
Wheel  width,  cm 
Cohesion,  N/cm 

3 

Penetration  resistance  gradient,  N/cm 
Instantaneous  center  of  rotation 
Displacement  or  distortion,  in.* 

Horizontal  soil  deformation  module,  cm  or  in. 

Cohesive  modulus  of  sirkage,  Ib/in. ^  ^ 

Friction  modulus  of  sinkage,  Ib/in. 

Lengths  of  lunule  zone,  bulldozing  zone,  and  zone  of  lateral 
flow,  respectively,  cm 

Factor  of  proportionality,  dimensionless 

Torque,  cm-N 

Sinkage  exponent,  dimensionless 
Wheel  radius,  cm 
Slip,  percent 

Shear- to- normal  stress  ratio,  dimensionless 
Mean  shear- to- normal  stress  ratio,  dimensionless 
Sliding  velocity,  cm/sec 

Widths  of  lunule  zone,  bulldozing  zone,  and  zone  of  lateral 
flow,  respectively,  cm 

Wheel  load,  N 

Distance  from  the  leading  edge  to  a  considered  point  on  the 
wheel- soil  interface,  in. 


*  The  dimensions  inches  and  pounds  refer  to  formulas  quoted  from  Anglo- 
Saxon  literature  (see  Conversion  Factor  table  on  page  ix).  All  expres¬ 
sions  pertaining  properly  to  this  report  are  in  the  metric  system. 


y  Coordinates  across  the  wheel-soil  interface,  measured  from 
the  center  line,  cm 

z  Sinkage,  cm 

Ot  Angle  at  the  wheel  center  between  the  leading  edge  of  the 
wheel-soil  interface  and  the  vertical  reference  line,  deg 
or  radians 

Oj  Angle  designating  the  separation  point  between  bulldozing 
and  trailing  zone,  according  to  Sela's  theory,  deg 

P  Angle  at  the  wheel  center  between  the  trailing  edge  of  the 
wheel- soil  interface  and  the  vertical  reference  line,  deg 
or  radians 

7  Inclination  of  the  shear  stress  vector  against  the  vertical 
plane  of  travel,  deg 

5  Position  angle,  angular  coordinate,  radians 

X  Experimental  constant,  dimensionless 

2 

0  Tangential  or  shear  stress,  N/cm 

■t  Normal  or  radial  stress,  N/cm^ 

<f>  Angle  of  internal  friction,  deg 
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CONVERSION  FACTORS,  METRIC  TO  BRITISH  AND  BRITISH  TO 
METRIC  UNITS  OF  MEASUREMENT 


Metric  units  of  measurement  used  in  this  report  can  be  converted  to  British 
units  as  follows: 


Multiply 

By 

To  Obtain 

centimeters 

0.3937 

inches 

meters 

3.2808 

feet 

Newtons 

0.2250 

pounds 

Newtons  per  square  centimeter 

1.4503 

pounds  per  square  inch 

penetration  resistance  gradient 
in  Newtons  per  cubic  centimeter 

11.0253 

0-6  in.  cone  index 

British  units  of  measurement  used  i.i  this  report  can  be  converted  to  metric 


units  as  follows: 

Multiply 

By 

To  Obtain 

inches 

2.54 

centimeters 

pounds 

4.444 

Newtons 

pounds  per  square  inch 

0.689476 

Newtons  per  square  centimeter 

SUMMARY 


"horrent  theoretical  concepts  of  the  soil-running  gear  interaction 
are  based  on  empirical  pressure- sinkage  and  shear  stress-displacement  re¬ 
lations.  It  was  found  that  the  tests  from  which  these  relations  were 
obtained  present  a  poor  analogy,  at  best,  to  the  soil-running  gear  inter¬ 
action.  In  particular,  the  use  of  the  shear  stress- displacement  relation 
for  an  analytical  evaluation  of  the  shear  stresses  at  the  contact  surface 
of  rigid  wheels  on  sand  was  found  to  be  misleading,  primarily  because 
equations  obtained  by  the  bevameter  or  dragged  plate  tests  describe  the 
soil  behavior  in  the  intermediate  state,  which  is  a  state  between  the 
elastic  and  plastic  states.  Recent  publications  point  out  that  part  of 
the  soil  in  the  immediate  vicinity  of  a  powered  wheel  is  in  a  state  of 
actual  plastic  flow  and  part  is  in  a  quasielastic  state.  Thus,  the  soil- 
bevameter  interaction  is  of  a  fundamentally  different  nature  than  the 
soil-wheel  interaction  which  makes  results  of  such  tedts  unsuitable  for 
predicting  wheel  performance. 

Sela's  theory  of  the  relation  between  a  rigid  wheel  and  dry  sand  is 
based  exclusively  on  the  shear  stress-displacement  concept,  and  thus  pro¬ 
vided  an  excellent  means  of  checking  the  concept  as  a  whole.  Using  a 
simple  approximate  relation  between  the  H/ffl  ratio  and  the  mean  shear- 
to-normal  stress  ratio  (t/<j)  taken  over  the  total  contact  surface,  the 
theory  was  checked  by  means  of  constant- slip  rigid  wheel  tests.  The 
experimental  findings  strongly  supported  the  conclusion  that  shear  stress- 
displacement  relations  are  irrelevant  to  soil-wheel  mechanics. 

A  new  theory  was  developed  that  attempts  to  assess  the  variation  of 
t/ct  along  the  soil-wheel  interface,  without  referring  to  shear  stress- 
displacement  relations.  The  proposed  theory  uses  recent  investigations  of 
the  soil  rupture  pattern  beneath  rigid  wheels  on  cohesionless  material  as 
the  basis  for  subdividing  the  total  soil-wheel  interface  into  three  zones 
of  different  soil  behavior.  Simple  assumptions  then  are  made  of  the  soil- 
to- wheel  movements  within  these  zones  to  evaluate  t/cj  .  The  approach 
developed  herein  was  checked  against  experimental  results  and  was  in  much 
better  agreement  with  the  test  data  than  were  previously  developed  theories. 
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CONTRIBUTION  TO  THE  MECHANICS  OF  RIGID  WHEELS  ON  SAND 


PART  I:  INTRODUCTION 


Background 


1.  The  forces  that  move  an  earthbound  vehicle  over  the  ground  Are, 
in  the  last  analysis,  soil  reaction  forces.  The  engine  serves  only  to 
generate  these  soil  reactions  by  transmitting  a  certain  mechanical  energy 
to  the  running  gear.  Recognition  of  this  fact  makes  the  interaction  be¬ 
tween  the  soil  and  the  running  gear  the  key  problem  of  theoretical 
vehicle- soil  mechanics. 

2.  Current  theoretical  approaches  to  this  fundamental  problem  are 
based  essentially  on  the  idea  of  stress-displacement  relations  for  soils. 
The  most  commonly  used  are  the  pressure- sinkage  relation  expressed  by  the 
equation 


a 


(1) 


and  the  shear  stress-displacement  relation 


t  =  (c  +  a  tan  0)(l  -  e"*^)  =  x  (l  -  e~^)  (2) 

luoJt 

Both  expressions  have  been  proposed  by  Bekker,1  the  latter  in  a  somewhat 

more  complex  form  to  account  for  peak  shear  stresses.  Equation  2  as  given 

.  .  P 

nere  is  a  simplification  introduced  by  Janosi  and  Hanamoto. 

3.  Although  it  is  obvious  from  the  manner  of  derivation  that  these 
equations  are  simply  mathematical  formulations  of  experimental  curves  ob¬ 
tained  under  special  conditions,  certain  researchers  have  tended  to  con¬ 
sider  them  as  laws  of  the  material.  Reference  to  these  equations  as 
"stress-strain  relations  for  soils"  or  to  some  of  the  parameters  as  "soil 
constants"  has  encouraged  this  tendency. 

4.  The  relations  stated  in  equations  1  and  2,  or  similar  ones,  are 
not  used  in  soil  mechanics,  despite  the  fact  that  knowledge  cf  the 


stress-strain  characteristics  for  soils  is  highly  desirable  for  many  soil 
mechanics  problems.  Because  the  problems  involved  are  important,  and  be¬ 
cause  the  concept  is  new,  stress- strain  relations  for  soils  should  have 
been  considered  with  studious  caution.  However,  it  appears  that  these  re¬ 
lations,  or  their  basic  ideas,  were  accepted  immediately  by  most  research 
workers  in  the  field  of  terrain- vehicle  mechanics.  Most  of  the  literature 
on  this  subject  is  related  to  secondary  issues  such  as  curve-fitting  tech¬ 
niques,  form  and  size  of  test  equipment,  establishment  of  a  definite  set 
of  "soil  constants,"  and  more  or  less  far-reaching  modifications  of  the 
analytical  expressions.  The  main  i.sue — whether  these  relations  are  appli¬ 
cable  at  all  in  the  field  of  soil-vehicle  mechanics — has  rarely  been  ques¬ 
tioned  seriously. 

Purpose 

5.  The  purpose  of  this  study  was  threefold: 

a.  To  critically  examine  the  shear  stress- displacement  rela¬ 
tions  for  soils  and  determine  whether  these  relations  are 
applicable  to  wheels  on  cohesionless  soils. 

b.  To  develop  new  means  of  experimentally  examining  any  theory 
on  the  ratio  of  shear-to-normal  stress  at  the  wheel- soil 
interface,  since  present  methods  are  not  very  reliable. 

c.  To  develop  a  theoretical  concept  regarding  the  shear-to- 
normal  stress  ratio  at  the  wheel- soil  interface. 

In  a  more  general  sense,  the  purpose  of  this  study  was  to  stimulate  a 
critical  reexamination  of  the  broadly  accepted  stress-displacement 
relations. 


Scope 

6.  Both  a  theoretical  and  an  experimental  approach  were  used  in 
this  study  to  investigate  the  distribution  of  the  tangential-to-normal 
stress  ratio  along  the  wheel-scil  interface  of  rigid  wheels  in  sand.  The 
experimental  portion  included  constant- slip  tests  in  sand  using  two  rigid 
wheels  instrumented  to  measure  tangential  and  normal  stresses  at  the 


2 


wheel- soil  interface,  and  the  length  of  the  interface  (contact  length). 
Wheel  load  and  width,  soil  strength,  and  slip  were  varied  to  determine 
their  effects  on  magnitudes  and  variations  of  the  stresses,  as  well  as 
on  the  performance  of  the  wheel.  Emphasis  was  placed  on  correctly  deter¬ 
mining  tne  M/RW  ratio  as  a  function  of  contact  length.  The  influence 
of  the  shear  stress  and  normal  stress  distribution  on  the  characteristic 
M/RW  ratio  was  evaluated  with  a  calculation  performed  by  computer. 
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PART  II:  THEORETICAL  DEVELOPMENT 


f-'  V" 


Review  of  Theories 

?.  It  must  be  emphasized  that  specific  stress- strain  relations  for 

■a 

soil  are  not  yet  established.  It  was  not  until  very  recently  that  first 
steps  were  taken  by  soil  mechanics  researchers  to  develop  theoretically 
founded  stress- strain  relations  for  soils,  in  the  sense  of  a  constitutive 
relation  of  the  material,  because  of  the  formidable  difficulties  inherent 
in  such  investigations. 

8.  In  soil  mechanics,  tne  soil  is  generally  idealized  and  assumed 
to  behave  like  an  elastic  or  plastic  material  or  according  to  some  other 
constitutive  relation,  such  as  viscoelasticity.  It  has  been  claimed  that 
this  approach  is  inadequate  in  soil- vehicle  mechanics;  that  a  "new  approach 
is  needed  to  relate  vertical  soil  deformation  to  ground  pressure  and  hori¬ 
zontal  soil  deformation  to  shear  strength."^  It  will  be  shown  later  that 
this  new  approach  aims  primarily  at  a  stress  domain  where  the  soil  is 
somewhere  between  the  elastic  and  the  plastic  state. 

9.  Another  "in  between"  aspect  of  the  problem  was  pointed  out  by 

5 

Deresiewicz,  who  attributed  to  granular  aggregates  a  stfLte  of  matter  be¬ 
tween  the  solid  and  the  liquid  state.  If  this  is  correct,  it  certainly 
presents  additiona.  difficulties  to  the  solution. 

10.  The  development  of  specific  stress- strain  relations  for  soils 
seems  possible  only  from  a  thorough  theoretical,  approach  that  takes  into 
account  the  discrete,  noncontinuous  character  of  the  material.  Such  a 
"physics  of  granular  aggregates"  has  recent ly  been  claimed  by  Kezdi,^  and 
basic  research  along  this  line  already  has  been  started. 

11.  Besides  Deresiewicz' s  review  of  previous  work  on  this  matter, 
there  are  other  recent  contributions  generally  characterized  by  the  attempt 
to  take  the  necessary  step  from  regular  packings  of  spheres  to  random  ar¬ 
rangements  of  odd-shaped  particles.  Green  and  Rivlin?  attempted  to  relate 
the  mechanics  of  structured  materials  to  that  of  continuous  media  by  intro¬ 
ducing  the  notion  of  multipolar  vector  fields.  Litwiniszyn  obtained  a 
fairly  accurate  description  of  sand  deformations  under  certain  conditions 


k 


of  flow  by  considering  the  discrete  movement  of  particles  as  a  stochastic 
process.  ( Berry®  gives  a  discussion  of  and  complete  reference  for 
Litwiniszyn' s  work.)  Murayama9  considered  the  probability  of  finite 
grain-to-grain  sliding  in  the  so-called  elastic  state  and  thus  explained 
the  relatively  low  values  of  Young's  modulus  for  sands.  Starting  from 
simple  assumptions  based  on  the  statistical  force  distribution  for  par¬ 
ticulate  matter,  Smoltczyk^®  developed  a  theoretical  approach  to  the  prob¬ 
lem  of  earth  pressures  prior  to  the  state  of  failure  (intermediate  state) 
that  may  become  a  link  between  the  given  elastic  state  at  rest  and  the 
state  of  plastic  failure.  Wiendieck^  investigated  the  structure  of  gran¬ 
ular  media  from  a  statistical  point  of  view. 

12.  These  investigations  demonstrate  the  multiple  problems  that 
have  yet  to  bj  overcome.  Once  a  thorough  understanding  of  the  internal 
deformation  mechanism  for  irregular  aggregates  is  gained,  perhaps  stress- 
strain  relations  can  be  established  by  integrating  the  elementary  process 
over  a  certain  volume.  Meantime,  the  quoted  stress-displacement  relations 
(equations  1  and  2)  must  be  taken  for  what  they  are:  analytically  ex¬ 
pressed  test  results,  with  no  theoretical  foundation. 

Applicability  of  the  Shear  Stress-Displacement 
Relations  to  Wheels  in  Sand 

4 

13.  In  a  comparative  study.  Wills  found  that  the  general  form  of 
the  shear  stress-displacement  curve  for  sand  obtained  with  a  torsional, 
ring  shear  apparatus  differed  significantly  from  that  obtained  with  a 
linear  shear  apparatus,  but  that  the  form  of  the  curve  is  not  dependent 
on  size  for  either  device.  He  further  noted  that  the  methods  (he  mentions 
three)  used  to  analyze  these  curves  are  highly  subjective  and  do  not  al¬ 
ways  yield  comparable  results.  For  instance,  in  a  given  case,  the  ratios 
of  corresponding  "horizontal  soil  deformation  moduli  k"  varied  from  1:2 
to  2:1,  and  the  absolute  magnitude  of  the  modulus  obtained  for  a  given  de¬ 
vice  depended  strongly  on  the  curve-fitting  technique  applied.  The  values 
are  as  follows: 

I 
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Method  of 
Analysis 

Ring  Shear 
Apparatus,  k 

St 

Linear  Shear 
Apparatus,  kg 

Ratio 

Bekker 

!-l/2 

3A 

1:2 

J&nosi 

1/2 

2:1 

Adams 

1-1/2 

3A 

1:2 

Equally  contradictory  results  are  obtained  when  ring  shear  and  linear  test 
results  are  compared  using  the  same  curve- fitting  techniques.  This  con¬ 
tradiction  is  remarkable  in  that  the  results  are  from  pure,  although 
slightly  different,  sliding  processes. 

14.  In  view  of  these  data,  the  application  of  shear  stress- 
deformation  relations  to  wheels  can  hardly  be  defended,  since  the  action 
of  a  rolling  wheel  on  the  soil  is  far  from  being  a  pure  sliding.  Their 
inapplicability  is  emphasized  by  closely  examining  the  physical  process 
that  occurs  in  the  sand  under  a  wheel  and  a  sliding  device,  respectively. 

15.  A  typical  shear  stress-displacement  curve  obtained  from  a 
horizontally  dragged  plate  test  is  shown  in  fig.  1.*  During  such  a  test, 

r  (psi)  the  soil  apparently  is  transformed 

from  a  presumably  elastic  state 
at  rest  into  a  perfect  plastic 
condition,  once  the  horizontal  as¬ 
ymptote  is  reached.  The  aim  of 
the  test  was  to  explore  not  the 
elastic  or  plastic  states,  but  the 
undefined  state  between,  which 
will  be  called  the  intermediate 

Fig.  1.  Shear  stress-displacement  re-  state  in  this  report.  This  process 

lation  obtained  with  linear  apparatus12  _  ,  „  .  .  . . 

of  transformation  necessarily  is 

accompanied  by  the  progressive  forming  of  a  certain  pattern  of  rupture  sur¬ 
faces  within  the  soil.  There  is  no  pattern  for  the  elastic  condition; 
there  is  a  completely  formed  one  for  the  plastic  condition. 

16.  What  does  a  partially  developed  rupture  pattern  look  like?  It 
might  be  called  an  incipient  rupture  pattern,  which  has  some  resemblance 

*  A  table  for  converting  British  units  (such  as  those  used  in  fig.  l)  to 
metric  units  is  given  on  page  ix. 
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to  the  definite  pattern  of  the  plastic  state  into  which  it  develops.  The 
main  feature  of  both  the  definite  and  incipient  rupture  patterns  is  a 
horizontal  critical  surface  somewhere  beneath  the  dragged  plate,  probably 
immediately  under  the  plate  or  along  its  grouser  tips,  when  the  plate  is 
provided  with  them.  Reece^  assumed  two  inclined  rupture  surfaces  that, 
together  with  the  plate,  confine  a  triangular  soil  wedge.  This  assumption 
has  no  bearing  on  the  following  discussion,  in  which  a  horizontal  rupture 
surface  was  assumed  for  the  sake  of  simplicity.  The  notion  of  stresses, 
as  applied  in  equation  2,  certainly  refers  to  a  horizontal  rupture  surface. 

17.  When  a  sliding  device  is  applied  to  the  soil,  shear  stresses 
at  the  incipient  horizontal  rupture  surface  build  up  gradually  from  zero 
to  Tmax  as  the  soil  progresses  from  the  elastic  condition,  through  an 
•undefined  intermediate  state,  to  the  plastic  condition.  For  the  shear 
stress-displacement  relation  to  apply  to  wheels,  it  is  necessary,  there¬ 
fore,  that: 

a.  The  wheel- soil  interface  be  an  incipient  rupture  surface. 

b.  The  soil  in  the  immediate  vicinity  of  the  interface  be  in 
neither  a  plastic  nor  an  elastic  state,  but  in  an  inter¬ 
mediate  one.* 
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18.  Boucherie  showed,  in  recent  investigations,  that  precisely 
the  contrary  happened  under  a  rigid  wheel  moving  on  an  artificial  two- 
dimensional  granular  medium  formed  by  a  random  packaging  of  aluminum  rolls 
(Schneebeli  model).  A  rigid  wheel  rolling  at  approximately  13  percent 
slip  is  shown  in  fig.  2,  and  a  schematic  diagram  of  the  phenomenon  in 
fig.  3.  At  the  forward  part  of  the  interface,  the  bow  wave  can  be  recog¬ 
nized.  The  soil  within  the  bow  wave  is  certainly  in  a  perfect  plastic 
condition.  The  trajectories  of  the  soil  particles  at  the  interface  are 
almost  perpendicular  to  it,  thus  the  rolling  surface  of  the  wheel  in  this 
region  is  assumed  to  be  a  real  rupture  surface. 

19.  At  the  rearward  end  of  the  interface,  a  moon-shaped  part  of  the 


*  Perfect  plastic  zones  and  real  rupture  lines  could  only  appear  at  the 
end  of  a  long,  intermediate  process,  i.e.  when  the  contact  surface  is 
long  enough  to  permit  the  soil  to  reach  plasticity.  The  statements  hold, 
therefore,  at  least  for  the  main  part  of  the  contact  zone. 


7 


Fig.  2.  Rupture  pattern  beneath  rigid  wheels^ 


DIRECTION  OF  TRAVEL 


Fig.  3.  Schema  of  rupture  pattern  beneath  rigid  wheels 
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soil,  which  has  been  termed  "lunule , "  seems  to  be  solidly  attached  to  the 
wheel  and  to  turn  around  the  instantaneous  center  of  rotation  together  with 
the  wheel.  Although  the  individual  soil  particles  in  the  formation  of  the 
lunule  change  constantly  as  the  wheel  moves  forward,  it  seems  that  the 
instantaneous  existing  lunule  behaves  like  a  solid  body.  Hence,  the  soil 
along  the  contact  zone  within  the  domain  of  the  lunule  has  to  be  considered 
in  an  elastic  state,  which  means  that  the  wheel- soil  interface  within  the 
domain  of  the  lunule  is  not  a  rupture  surface. 

20.  This  general  schema  is  valid  through  the  whole  positive  slip 
range;  only  the  relative  proportions  of  the  two  regions  change.  At  zero 
slip,  the  forward  bulldozing  zone  is  predominant;  at  100  percent  slip, 
the  outer  limit  of  the  lunule  coincides  with  the  wheel  circumference,  thus 
covering  the  total  contact  length.  Although  Boucherie's  experimental 
study  was  only  two-dimensional,  there  is  no  reason  to  doubt  that  the  same 
general  pattern  of  soil  kinematics  exists  under  an  actual,  wheel.  This  is 
confirmed  by  recent  research  by  Wong  and  Reece, ^  who  made  photographic 
studies  similar  to  Boucherie's  of  the  coil  movements  beneath  a  wheel,  but 
used  an  actual  sand  fill,  rather  than  an  artificial  medium.  These  inves¬ 
tigations  showed  that  the  pattern  under  an  actual  wheel  is  the  same  as 
the  one  shown  in  figs,  2  and  3* 

21.  Since  no  proof  has  been  given  that  the  lunule  behaves  like  an 
elastic  solid,  it  might  be  argued  that  the  soil  in  this  region  is  in  the 
intermediate  state.  However,  this  would  require  that  the  intermediate 
state  develop  behind  the  plastic  state  (measured  from  the  leading  edge  of 
the  wheel-soil  interface)  which  contradicts  the  way  shear  stress- 
displacement  relations  are  usually  applied. 

22.  The  existence  of  quasisolid  soil  masses  attached  to  rigid  ele¬ 
ments  moving  into  or  through  cohesionless  soils  is  well  known  in  soil 
mechanics,  Biarez^  observed  such  rigid,  nondeforming  soil  bodies  in 
almost  all  of  his  investigations  of  various  models  of  construction  elements 
in  the  plastic  soil  condition  and  demonstrated  it  using  a  special  photo¬ 
graphic  technique  that  allowed  the  camera  to  move  with  the  elemertt.  The 
general  shape  of  the  soil  bodies  was  determined  by  the  conditions  of 
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ACTIVE  PLASTIC  ZON 


PTURE  LINE  OF  THE 
PILE-BALL  SYSTEM 


kinematic  compatibility,  while  their  actual  size  seemed  to  be  such  that  the 
effort  needed  to  move  the  elements  was  a  minimum. 

23.  A  pivoting  pile  bears  a  striking  resemblance  to  the  wheel  prob¬ 
lem.  As  shown  in  fig.  4,  a  ball-socket  kind  of  hinge  develops  in  the  lower 

part  of  the 

. .  ^  (APPLIED  FORCE) 

«  medium  and  has 

/ 

^  the  same  function 

/  as  the  lunule, 

A  i.e.,  to  make  its 

rotation  kine¬ 
matically  possi¬ 
ble  and/or  to 
facilitate  the 

passive  plastic  zone  rotation  around 

the  instanta¬ 
neous  center  of 

pile-ball  system  rotation.  P  .as- 

tic  soil  masses 
appear  only  in 
the  upper  part  of 
the  medium.  This 
schema,  in  com¬ 
bination  with  the 

principle  of  min- 

Fig.  4.  Pivoting  pile  showing  ball-socket  hinge  .  „ _ „ 

imum  energy  or 

rotation,  was  used  in  a  new  approach'  to  the  stability  of  sheet  pile  walls 
and  matched  the  results  of  model  tests  better  than  previous  theories. 

24.  A  similar  approach  to  the  mechanics  of  a  rolling  wheel  might 
clarify  some  open  questions.  But  before  this  can  be  undertaken,  some 
light  must  be  shed  on  the  stress  distribution  along  the  line  ABC  in  fig.  3. 
The  apparently  rigid  lunule  might  then  appear  as  a  key  feature  of  the 
kinematics  of  wheel  movement. 

25.  From  the  foregoing  (see  paragraph  17),  it  can  be  stated  that: 

a.  A  rupture  surface  occurring  at  the  wheel-soil  interface  is 
”  fully  developed. 


Fig.  '*■.  Pivoting  pile  showing  ball-socket  hinge 
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b.  The  soil  along  the  interface  beneath  a  rigid  wheel  is 

either  in  a  state  of  perfect  plastic  flow  (bulldozing)  or 
in  a  quasielastic  state  (lunule). 

Clearly,  these  conditions  do  not  meet  the  theoretical  requirements  implied 
by  the  use  of  equation  2. 

2 6.  It  has  been  demonstrated  that  the  action  of  a  wheel  upon  co- 
hesionless  soil  is  completely  different  from  that  of  sliding  devices.  The 
shjar  stress-displacement  equation,  therefore,  cannot  be  applied  to  wheels 
in  sand. 

Brief  Outline  of  Sela's  Theory 
17 

27.  Sela's  theory  is  based  entirely  on  the  concept  of  shear  stress- 
displacement  relations,  and,  as  such,  provides  an  excellent  opportunity  to 
check  the  basic  theoretical  concept. 

28.  For  cohesionless  soils  (c  =  0),  equation  2  may  be  rewritten: 

t  =  |  =  tan  0(1  -  e‘J^)  (3) 

Usually  applied  to  wheels  or  tracks 

j  =  xs  (4) 

based  on  the  assumption  that  "the  distortion  j  might  be  expressed  as  a 
product  of  slip  (s)  and  the  distance  (x)  measured  from  the  beginning  of  the 
contact  area  to  Ihe  place  where  the  shearing  strength  of  the  soil  is  devel¬ 
oped."  Expressed  in  radians  (see  fig.  8,  page  17) 

x  =  R(a  -  8)  (5) 

Accordingly: 

t  »  i  =  tan  <j>  [l  -  e'R^a  '  5^s/k]  (6) 

17 

29.  Sela  refined  equation  6  to  take  into  account  the  effect  of  the 
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bow  wave,  but  neglected  lateral  flow.  For  the  bow  wave  zone  he  proposed 
the  equation 


t  =  I  =  tan  f  jl  -  e-R[(a  "  6>s  +  '  “dWAl  (7) 

where  the  angle  designates  the  separation  point  between  the  bow  wave 
region  (as  defined  by  Sela)  and  the  rear  part  of  the  interface,  and  m  is 
a  factor  of  proportionality.  The  angle  is  defined  by  the  relation 

tan  ad  =  ^  =  tan  $  [l  -  e"1^  '  5^s/k] 


30.  An  example  of  the  resulting  theoretical  t/ct  curve,  together 

with  the  experimental  measurements  made  with  pressure  cells,  is  reproduced 
17 

from  Sela’s  work  in  fig.  5*  According  to  Sela's  theory,  the  variation 
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Fig.  5.  t/ct  measurements 
using  Sela’s  theory 


of  the  t/c  ratio  over  the  contact  length  is  generally  characterized  by 
an  increase  toward  both  ends  of  the  wheel- soil  interface,  with  the  lowest 
value  at  the  point  of  separation  between  the  forward  bulldozing  and  the 
rearward  trailing  zone. 


Experimental  Ways  to  Check  Theoretical  Concepts 
of  Shear /Normal  Stress  Ratio 


Tangential  and  normal 
pressure  cell  recordings 

31.  To  use  this  method,  a  set  of  carefully  calibrated  and  protected 
pressure  cells,  placed  flush  with  the  rolling  surface  of  the  wheel,  re¬ 
cords  the  tangential  and  normal  stresses  as  the  wheel  surface  moves  through 
the  soil.  Sela^f  used  this  method  to  check  his  theoretical  concept  and 
found  that  the  general  shapes  of  the  distribution  curves  for  the  measured 
and  predicted  variation  of  the  t/<~  ratio  matched  reasonably  well. 

32.  However,  pressure  cell  recordings  are  not  completely  reliable, 
since  neither  the  influence  of  the  protective  membrane  necessary  to  pre¬ 
vent  sand  wedging  against  the  cell  nor  the  possible  arching  effects  are 
known.  In  WES  tests  using  cell  recordings,  discontinuous  rupture  phenomena 
in  the  sand  caused  considerable  vibration  of  the  wheel  that  made  the  re¬ 
corded,  stress  curves  completely  unreadable  in  some  instances.  The  tangent 
tial  stress  records  were  affected  more  severely,  and  the  degree  of  dis¬ 
turbance  seemed  to  increase  as  slip,  soil  strength,  and  load  increased. 

When  actual  t  and  a  readings  were  used  to  determine  the  pertinent  t /a 
ratios,  an  enormous  scatter  resulted  (see  figs.  6  and  7*).  The  data  scat¬ 
ter  can  be  reduced  either  by  smoothing  the  original  curves  by  hand  or  by 
electronically  filtering  the  signals  from  the  pressure  cells.  The  first 
method  introduces  subjective  personal  judgment;  the  second  usually 

*  In  these  figures,  some  of  the  most  scattered  points  are  not  shown.  To 
obtain  some  representable  results  in  the  high  slip  range,  it  was  neces¬ 
sary  to  smooth  the  t  and  a  recordings  by  hand.  These  data  represent 
the  stresses  measured  directly  beneath  the  axle.  When  the  ratio  is  com¬ 
puted  near  the  midpoint  of  the  wheel- soil  interface,  the  scatter  is  con¬ 
siderably  less. 
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entails  a  certain  time  delay  unless  highly  sophisticated  equipment  is 
used.  Both  methods  tend  to  overlook  the  fact  that  these  disturbing  vibra¬ 
tions  are  a  real  phenomenon  and  not  just  a  nuisance  to  be  eliminated. 

33*  Also  unreliable  is  the  primitive  method  of  using  much  less 
sensitive  pressure  cells  to  circumvent  these  problems  (apparently  the  case 
in  Sela's  measurements),  especial ly  at  the  extremes  of  the  wheel- soil 
interface,  where  "readings  of  the  shear  and  normal  stresses  are  low  and 
thus  a  small  error  in  reading  causes  a  large  error  in  the  x /a  value. 
Because  of  the  shortcomings,  pressure  cells  give  only  a  rough  idea  of  the 
variations  and  the  order  of  magnitude  of  the  t/c  ratios. 


Comparison  of  mean  shear/ 
normal  stress  ratio  to  m/kw 

34.  This  indirect  method  of  checking  the  theoretical  concept  is 
based  on  the  consideration  of  the  mean  value  of  t/<j  (=1^)  taken  over  the 
whole  contact  area,  rather  than  on  its  variation  along  the  contact  length. 

It  can  be  shown  that  t^  is  roughly  equal  to  the  M/} RW  ratio.  Since  the 
tjj,  value  can  be  obtained  readily  from  any  x/a  theory  (by  graphical 
construction,  if  necessary),  and  can  be  easily  and  accurately  meas¬ 

ured,  the  relation  between  the  two  constitutes  a  convenient  tool  of  control. 

35-  The  ratio  can  be  expressed  by 


M_ 

KW 


where 


2 

r  a  b/2 

If  ' 

-p  0 

‘(B.y)  °(6,y)  “  ^ 

a  b 

fj 

-P  C 

l2 

ff(8}y)  cos  5  d6  dy 

+  2 

a  b/2 

hs.y)  °(6,y) sin  6  dB  * 

=;(8) 


o,„  \  =  normal  pressure  as  a  function  of  angle  6  and  the  distance 

'  y  from  the  center  of  the  rigid  wheel  (fig.  8) 

^(s  y)  =  corresPondins  T/h  f'uhcti011 

In  equation  8,  the  second  term  of  the  denominator,  which  represents  the 
contribution  of  shear  stresses  to  the  load-carrying  capacity  of  the  wheel, 
is  a  very  small  quantity  compared  with  the  first  term,  which  represents 
the  contribution  of  normal  stresses,  a  .  Thus,  if  the  second  term  is  ne¬ 
glected,  only  a  small  negative  error  is  generated,  and  this  can  be  roughly 
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Fig.  8.  Rigid  wheel  coordinates 

compensated  by  simultaneously  cancelling  the  cosine  factor  in  the  first 
term,  which  slightly  increases  the  value  of  the  denominator. 

36.  For  a  qualitative  evaluation,  equation  8  therefore  can  be  sim¬ 
plified  to  1-  -1 

a  b/2 


M  2  V.yj'W)®* 

P  a  b/2  1 

2\^{  a(E,y)  d6 


Equation  9  can  be  simplified  further  by  replacing  the  function 
with  its  mean  value  mathematically  defined  by 

f  a  b/2  "I 


14  4  Wd5dy 


m  b(a  +  P) 

37.  Since  tm  is  constant  for  a  given  contact  zone,  it  can  be  put 
before  the  integral  sign,  and  the  two  remaining  integrals  in  the  numerator 
and  denominator  cancel.  This  leads  to  the  simple  approximation  formula: 


*L«t 

KW  m 


This  constitutes  a  rather  crude  evaluation,  of  course,  but  the  general 
validity  of  the  approximate  formula  is  remarkable  since  it  can  be  devel¬ 
oped  without  considering  wheel  width,  slip,  or  soil  strength.  The  devi¬ 
ation  from  the  true  1^/RW  ratio  depends  only  on  the  t  and  a  functions 
and  on  the  integration  limits  of  a,  and  P  . 

38.  To  check  the  degree  of  approximation,  equation  11  was  compared 
with  the  exact  equation  8,  which  had  been  solved  by  a  computer  for  all 
possible  combinations  of  five  t  distributions  and  four  a  functions. 
Assuming  P  =  0.2  a  ,  which  is  close  to  its  experimental  value,  ty/HW  was 
calculated  for  a  values  varying  from  0.1  to  1.1  radians.  Variations  in 
the  cr  and  t  functions  over  the  width  of  the  wheel  can  be  neglected 

in  these  hypothetical  computations. 

39.  The  assumed  a  and  t  distributions  along  the  contact  length 
are  tabulated  in  fig.  9*  Although  they  are  all  distorted  sine  and  cosine 


t  -  distributions 


a  -  distributions 
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Fig.  9*  o  and  t  distributions  used  for  computations 


functions,  the  variety  of  their  main  features  represents,  in  a  qualitative 


18 


-  tit  ~j  " 


manner,  all  possible  t^es  of  distributions  hat  can  reasonably  be  assumed 
for  rigid  wheels  in  sand.  It  is  the  type  of  distributions  that  is  impor¬ 
tant,  not  how  they  are  mathematically  expressed.  The  t  functions  have 
been  so  formulated  that  their  tm  value,  as  defined  by  equation  10,  is 
equal  to  unity.  By  multiplying  them  by  a  constant  n  ,  any  desired  tm 
value  can  be  introduced  (n  •  1  =  i^).  The  results  of  these  computations 
are  listed  in  table  1,  and  some  are  graphically  illustrated  in  fig.  10. 


0.3  h  e 


0.825  t  -  0.2475 


0.1  U  n  . 


CONTACT  ANGLE 


Pig.  10.  Results  of  computations  for  selected  o-t  combinations 


These  were  chosen: 


a.  The  M/lN  curves  with  the  greatest  positive  or  negative 
~  deviation  from  the  value  (combinations  B2  or  C3  and 

S3»  respectively). 

b.  The  curve  for  the  simplest  c-t  combination  (Al). 

c.  The  curve  for  the  most  probable  o-t  combination  among 
those  investigated  (D5).* 

40.  Of  the  360  calculated  v/r W  values  (table  l),  only  10  deviated 
more  than  25  percent  from  tm  value.  These  axe: 


Combination  B2 


t  * 
m 

0.1 

a 

m/rw 

Deviation 

i 

0.7 

0.1258 

+25.8 

0.9 

0.1281 

+28.1 

1.1 

0.1312 

+31.2 

Combination  B3 


a 

Deviation 
H'HW  i 

0.5 

0.3702 

-26.0 

0.7 

0.3657 

-26.9 

0.9 

0.3634 

-27.3 

1.1 

0.3634 

-27.3 

Combination  C3 


t  = 
m 

0.1 

q 

h/rw 

Deviation 

* 

0.7 

0.1288 

+28.8 

0-9 

0.1354 

+35.4 

l.l 

0.1451 

+45.1 

However,  none  of  these  are  probable  combinations.  C3  is  the  most  improb¬ 


able  of  all  those  investigated,**  and  the  extreme  positive  deviations  occur 
only  for  high  contact  angles  (a  >  0.7)  at  a  very  low  tjj,  value.  But 
Sela's  theory,  as  well  as  the  new  concept  that  will  be  presented  (Part  V), 
indicates  that  for  high  contact  angles  the  value  is  greater  than  0.1, 
so  incidental  extreme  positive  deviations  can  be  disregarded.  The  highest 
negative  deviations  (B3)  are  barely  above  25  percent  and  can  also  be  dis¬ 


regarded.  Equation  11  can  then  be  written: 


As  mentioned  previously,  equation  11a  is  an  easy-to-handle  tool  of  control. 
If  the  experimental  m/RW  values  do  not  match  the  theoretical  tm  values 
within  a  25  percent  error  margin,  the  basic  theory  may  be  considered 
erroneous. 


*  The  t  function  (5)  comes  nearest  to  the  distribution  found  by  Sela. 
All  measurements  of  a  distributions  under  rigid  wheels  on  sand  indi¬ 
cate  a  type  of  distribution  comparable  to  D  (fig.  9). 

**  The  most  probable  combination,  D5,  presents  only  negative  deviations, 
with  the  extreme  value  of  -21.6  percent. 
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Qualitative  conparision  of 

variations  of  t  to  m/hW 
_ m _ Y 

4l.  If  it  is  assumed,  that,  during  a  particular  test  series,  the 
typical  features  of  the  t  function  and  the  a  distribution  do  not 
change,  equation  11  can  be  rewritten: 


KW 


A  t 
ot  m 


(12) 


The  subscripts  a  and  t  indicate  that  the  coefficient  A  is  valid  only 
for  one  particular  a-t  combination,  not  for  the  whole  group  of  possible 
combinations,  and  is  intended  to  approximate  only  one  individual  curve. 

For  example,  the  coefficient  A^  =  0.825  approximates  the  D5  combination 
+5$  (fig.  9)  if j  as  before,  the  extreme  deviations  at  a  =  0.9  (+7$)  and 
a  =  1.1  (+11.5$)  for  =  0.1  are  disregarded.  Since  a  greater  devia¬ 
tion  may  occur  in  the  actual  o-t  combination,  it  seems  fair  to  write: 

I? -  V.  ± 10*  ^ 

42.  The  effective  Agt  value  is  unknown  for  any  real  problem,  so 
aquation  12a  can  only  be  used  to  compare  the  variations  of  m/kW  and 
tm  ,  not  their  absolute  values.  Because  the  10  percent  error  margin  is 
sman,  equation  12a  can  be  stated: 


M_ 

KW 


t 

m 


(12b) 


This  means  that  M/KW  varies  directly  as  ,  if  it  can  be  assumed  rea¬ 
sonably  that  the  general  features  of  the  a  and  t  distributions  do  not 
change  radically  during  th*.  test  program. 


PART  III:  TEST  RESULTS  USED  TO  EVALUATE  THEORIES 

43.  A  series  of  constant-slip  tests,  in  which  the  contact  length 

was  varied  by  increasing  the  load,  provided  the  conditions  necessary 

(i.e.  o-t  combination  is  basically  invariable  with  contact  length)  for 

equation  12b  tc  be  applied.  These  tests  were  run  in  an  air-dry  sand, 

classified  SP  under  the  Unified  Soil  Classification  System.  Sand  at  two 

strength  levels  was  used,  and  the  penetration  resistance  gradients  were 

approximately  2.72  and  5.44  N/cm^  (0-6  in.  cone  indexes  of  30  and  60, 

respectively).  The  equipment  and  test  techniques  used  have  been  de- 

3.3 

scribed  by  McRae,  Powell,  and  Wismer. 

44.  The  normal  or  radial  stresses  were  measured  with  deflecting 
diaphragm- type  transducers  mounted  so  that  the  active  portion  of  each 
cell  was  flush  with  the  surface  of  the  center  line  of  the  wheel.  Tangen¬ 
tial  stresses  were  measured  with  a  cantilever  beam- type  transducer.  One 
end  of  the  beam  was  flush  with  the  outer  surface  of  the  wheel;  the  other 
securely  anchored  within  the  wheel. 

45.  The  contact  length  was  of  particular  importance  in  this  study. 
The  classical  measure  of  sinkage  was  not  adequate  for  determining  contact 
length,  since  it  does  not  account  for  the  bow  wave  before  and  the  up¬ 
heaval  behind  the  wheel.  Pressure  cells  were  considered  inadequate  to 
determine  the  contact  length,  because  of  the  difficulties  discussed  in 
paragraphs  32  and  33-  Therefore,  a  mechanical  switch  was  placed  at  the 
wheel  surface,  and  the  time- in- soil  recorded  on  a  direct-writing  oscillo¬ 
graph.  However,  a  force  of  approximately  3*5  N  was  needed  to  activate 
the  switch,  so  the  recorded  contact  length  was  somewhat  smaller  than  the 
actual  one,  but  since  it  did  not  affect  the  qualitative  interpretation  of 
the  test  results  (see  footnote  on  page  32),  no  effort  was  made  to  cor¬ 
rect  the  systematic  error. 

46.  m/hw  curves  as  functions  of  contact  length  are  plotted  for  four 
test  series  in  figs.  11- l4.  Since  slip,  soil  strength,  and  load  were  not 
the  same  as  the  preestablished  values  in  all  cases,  the  data  had  to  be  cor¬ 
rected  before  the  curves  could  be  drawn.  The  actual  measured  values  are 
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shown  by  marked  points;  the  areas  which  .r  ild  contain  the  corrected  read¬ 
ings  are  circumscribed  by  rectangles. 


CONTACT  LENGTH, CM 

Fig.  1...  Measured  ^/RW  versus  contact  length 
Mortar  sand 

G  =  5.44  N/cm3  tan  ^  =  0.65 
R  =  35.5  cm  b  =  30.5  cm 


Measured  ty/RW  versus  contact  length 
Mortar  sand 

=  5.44  N/cm3  tan  $  =  0.65 
=  35.5  cm  b  =  7.6  cm 
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PART  IV:  EVALUATION  OF  SELA’S  THEORY 


h'j.  Two  criteria  have  been  developed  for  comparing  ty/sw  curves  to 
the  theoretical  tm  curves.  The  first,  equation  11a,  states  that  for  the 
same  contact  length,  the  absolute  h/fW  values  should  coincide  with  the 
absolute  t^  values  within  a  25  percent  error  margin.  In  the  second  method 
of  checking,  using  equation  12b,  certain  qualitative  descriptive  terms 
should  apply  to  both  the  experimental  M/RW  curves  and  the  theoretical  tm 
curves. 

48.  The  main  qualitative  features  of  the  Vi/'B W  curves  in  figs.  11- 14 

are: 

a.  The  higher  the  slip,  the  higher  the  M/RW  ratio. 

b.  The  first  portion  of  each  w/hW  curve  is  concave  upward. 

c.  The  length  of  the  concave  upward  portion  of  the  M/RW  curves 
seems*  to  increase  as  the  slip  increases. 

d.  The  curvature  in  the  concave  upward  portions  of  the  m/rW 
curves  is  more  pronounced  at  lower  slip  values. 

e.  The  data  strongly  indicate  that  for  zero  contact  length, 
the  m/TsSU  ratio  is  not  zero,  but  a  definite,  positive 
value . 

Determination  of  the  tffl  Values  Using  Sela's  Theory 


49.  To  check  Sela's  theory  against  the  test  results,  the  theoretical 
tjjj  values  of  this  theory  must  be  determined.  Nine  different  t-distribntions 
were  obtained  by  applying  Sela's  theoretical  concept  to  the  WES  test  data. 
These  are  given  in  fig.  15  and  are  based  on  the  following  data: 

R  =  71/2  =  35.5  cm 
k  =  2.54  cm  (l  in. ) 
m  =  0.5 
tan  =  0.5 

slip  =  3 1o,  14$,  and  33$ 


*  Some  curves  were  too  short  or  too  flat  to  recognize  the  end  of  the  con¬ 
cave  upward  portion. 


k  and  m  are  not  measured  for  this  sand,  since  their  exact  determination 
is  highly  subjective. 5  The  values  used  by  Sela^  were  k  =  2.5k  cm  and 
m  =  O.65  ,  but  he  pointed  out  that  the  value  of  m  was  too  high;  there¬ 
fore,  a  slightly  smaller  value  (0.5)  was  used  in  this  analysis.  Tan  <f>  - 
0.5  was  chosen  as  a  convenient  reference  value.  The  actual  friction  co¬ 
efficients  of  the  test  sand  were  0.6  and  O.65,  respectively,  so  the  t 
and  t  values  (figs.  15  and  l6,  respectively)  were  multiplied  by  a  factor 


Fig.  16.  Variation  of  tm  with  contact  length, 
using  Sela's  theory 


of  1.2  and  1.3,  respectively.  The  pertinent  tffi  values  were  obtained  from 
individual  t-distribution  curves  by  graphical  construction  (fig.  15d)  and 
are  plouted  as  functions  of  contact  length  in  fig.  16. 


Theoretical  and  Experimental  Results  Compared 
Using  Sela's  Theory 

Comparison  of  absolute  values 

50.  The  corrected  tni  curves  of  fig.  lb  and  the  +25  percent  error 
limits  were  superimposed  onto  the  corresponding  I^/RW  curves  in  figs.  17 
and  18.  To  avoid  obscuring  the  graphical  representation,  the  I^/EW  curves 
were  regrouped  in  these  figures.  For  Sela's  theory  to  be  correct,  all 
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experimental  m/yN  curves  should  certainly  fall  within  the  +25  percent 
limits;  however;  to  a  significant  degree,  they  do  not.* 

Comparison  of  variations 

51.  The  second  check  of  the  theory,  provided  by  equation  12b,  states 
that  h/hW  varies  with  contact  length  in  the  same  manner  as  ,  if  the 
a- t  combination  is  assumed  to  be  basically  invariable  with  contact  length. 
This  assumption  is  essentially  valid  for  constant-slip  tests;  therefore, 
the  t^  curves  of  fig.  l6  should  match,  point  by  point,  the  main  qualita¬ 
tive  features  of  the  m/kW  curves  (paragraph  48)  if  Sela's  theory  is 
correct.  For  the  five  enumerated  points,  only  point  a  is  in  agreement. 
Sela's  theory,  therefore,  must  also  be  rejected  on  the  basis  of  the  second 
criterion. 


Discussion  of  Results  Obtained  by  Sela 

52.  In  Sela's  investigation,  the  measured  i/a  distribution  agreed 
closely  with  his  theoretical  predictions.  This  can  be  explained  by  com¬ 
paring  the  results  of  tests  on  the  weaker  sand  (tan  0  =  0.6,  fig.  17)  with 
those  on  stronger  sand  (tan  0  =  O.65,  fig.  18).  For  high  contact  angles 
(a  +  P  >  0.5  radians),  the  M03?  curves  for  the  weak,  sand  are  all  within 
the  limits  of  Sela's  theory;  for  the  dense  sand,  they  are  not.  Since  the 
sand  used  by  Sela  was  even  weaker  (tan  0  =  0.44,  bevameter  reading),  and 
the  contact  angles  (a  +  P)  ranged  from  0.6  to  1.2  radians,  it  is  not 
surprising  that  the  measured  t/ct  variations  matched  his  theory  fairly 
well.  The  lack  of  agreement  appears  only  when  a  broader  spectrum  of  con¬ 
tact  angles  and  soil  strengths  is  considered. 

Validity  of  Any  Shear  Stress  Relation 

53.  The  second  method  of  checking  disqualifies  all  theories  on  the 


*  This  percentage  would  not  change  significantly  if  the  error  due  to  the 
faulty  measurement  of  contact  length  were  eliminated.  The  displacement 
to  the  right  of  the  experimental  M./RW  corves  would  diminish  the  portion 
to  the  left  of  the  curves  and  increase  it  on  the  right. 


shear-to-normal  stress  ratio  of  wheels  developed  from  the  shear  stress- 
displacement  relation  (equation  2).  The  very  mathematical  nature  of  this 
exponential  function  will  not  permit  conditions  b  and  e  (paragraph  48) 
to  be  fulfilled  by  such  a  concept.  When  equation  2  is  used  as  a  basic 
relation  to  predict  the  i/o  ratio  for  wheels,  the  resulting  tju  curves 
as  functions  of  contact  length  will  always  be  concave  downward  and  will 
always  pass  through  the  zero  point  for  zero  contact  length,  unless  the 
equation  is  distorted  by  additional  terms  in  such  a  manner  that  its  ex¬ 
ponential  nature  can  no  longer  be  recognized. 
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PART  V:  PROPOSED  NEW  CONCEPT 


Development  of  Concept 


54.  It  has  been  shown  that  the  soil  rupture  patterns  beneath  moving 
wheels  and  those  obtained  by  simple  shear  test  devices  are  not  the  same. 
Therefore,  any  approach  to  soil-wheel  mechanics  must  be  based  entirely  upon 
consideration  of  an  actual  wheel  and  the  corresponding  pattern  of  soil 
flow.  The  rough  general  features  of  the  soil  deformation  pattern  beneath 
moving  rigid  wheels  have  been  established  through  research  based  on  visual 
observation  of  soil-wheel  interaction  conducted  at  Grenoble^  and  New¬ 
castle  upon  Tyne . ^  From  these  studies,  two  principal  regions  within  the 
wheel- soil  interface  have  been  distinguished:  a  rearward  region  in  con¬ 
tact  with  the  soil  in  the  elastic  state  (lunule),*  and  a  forward  region  in 
contact  with  the  soil  in  a  perfect  plastic  state  of  actual  flow. 

55-  These  studies,  however,  were  limited.  No  quantitative  data 
are  available  that  reflect  changes  in  the  lunule  due  to  slip,  wheel  geom¬ 
etry,  or  soil  properties.  The  shapes  of  the  three-dimensional  lunule  and 
forward  soil  body  are  unknown.  Accordingly,  the  principal  assumptions 
concern  the  boundaries  of  three  different  zones  within  the  wheel- soil  in¬ 
terface:  the  lunule  zone,  the  bulldozing  zone,  and  the  lateral  flow  zone, 

56.  In  the  three-dimensional  case,  it  is  indeed  logical  to  sub¬ 
divide  the  main  forward  region  into  two  zones  of  different  soil  behavior. 
The  principal  assumption  is  that  a  state  of  actual  soil  flow  is  maintained 
for  the  total  forward  region,  but  the  two  zones  differ  in  the  direction 


*  This  point  of  view  is  not  textually  supported  by  Wong  and  Reece. ^ 

While  Boucherie^  refers  to  the  moon- shaped  soil  mass  as  being  rigidly 
attached  to  the  wheel  and  moving  around  the  instantaneous  center  of  rota¬ 
tion,  Wong  and  Reece^  consider  it  as  being  in  an  active  (plastic)  state 
bounded  by  a  logarithmic  spiral.  Besides  the  misleading  term  "active 
state"  (an  active  state  is  defined  as  being  a  state  of  plastic  failure  in 
which  the  weight  of  the  considered  soil  mass  supports  the  failure  proc¬ 
ess,  which  is  not  the  case  here),  another  controversial  point  has  to  be 
mentioned.  Wong  and  Reece  noted  indeed  that  the  s:>il  "moves  backward  by 
means  of  rotation  approximately  about  the  instantaneous  center  of  move¬ 
ment,"  a  statement  which  is  somewhat  in  contradiction  to  the  hypothesis 
of  a  logarithmic  spiral  boundary. 
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of  the  soil  flow.  In  the  bulldozing  zone,  the  soil  particles  are  assumed 
to  move  forward  in  the  vertical  plane  of  travel;  in  the  lateral  flow  zone, 
the  soil  is  assumed  to  escape  from  under  the  wheel  by  roving  parallel  to 
the  vertical  plane  through  the  wheel  axle.  The  division  into  zones  of  the 
wheel- soil  interface  is  shown  in  fig.  19. 


Fig.  19.  Zoning  of  the  wheel-soil  interface 


Lunule  zone 

57*  The  width  of  the  lunule  zone  is  assumed  to  be  equal  to  the  width 
of  the  wheel: 


wL  =  b  (13) 

The  relative  length  of  the  lunule  zone,  with  respect  to  the  total  contact 
length,  is  primarily  a  function  of  slip,^  At  very  low  positive  slips,  the 
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lunule  is  barely  distinguishable ;  ^3  at  100  percent  slip,  it  covers  the 
whole  contact  length.^  It  is  assumed  that  the  relation  between  the 
length  of  the  lunule  zone  and  slip  is  linear.  The  portion  of  the  lunule 
zone  forward  of  the  vertical  reference  position  (5  =  0)  is  then  (fig.  19) 

R6*  =  sRa 

If  it  is  assumed  that  0  -  0.2a  ,  the  rearward  portion  of  the  lunule  zone 
has  to  be  limited  for  s  >  0.2  ,  which  is  approximately  the  maximum  ob¬ 
served  in  testing.  Hence 

RB"  =  0.2RQ:  (s  >  0.2) 

L 

For  slip  values  less  than  0.2,  it  is  assumed  that  RB"  =  R6^  .  The  total 

JU  ll 

length  of  the  lunule  zone  is  then 

iT  =  R(5*  +  b~)  =  2sRa  (s  <  0.2) 

1j  Jj  Jj 

iT  =  (0.2  +  s)Ra  (s  >  0.2)  (14) 

58.  The  x/a  ratio  (t)  for  that  section  of  the  wheel- soil  inter¬ 
face  within  the  lunule  is  less  than  tan  because  the  lunule  is  in  a 
nonplastic  state.  It  was  found  to  be  very  difficult  to  make  any  assump¬ 
tions  about  the  t- variation  within  this  zone,  even  though  the  t  value 
along  the  outer  limit  of  the  lunule,  which  is  a  rupture  surface,  is  known 
( t  =  tan  0).  Therefore,  the  t  value  was  measured  at  the  center  line 
of  the  wheel  by  means  of  normal  and  tangential  pressure  transducers  in¬ 
corporated  at  the  rolling  surface.  The  readings  for  5=0  are  given  in 
figs.  6  and  7;  the  large  scatter  inherent  in  sucn  measures  is  apparent. 

At  high  slip  values,  only  those  t  values  that  could  be  read  clearly  on 
the  recordings  are  listed.  Nevertheless,  a  general  trend  of  the  t^^ 
values  as  a  function  of  slip  can  be  recognized  that  can  be  expressed  as 
the  hyperbolic  function 

*  Because  of  the  sand-glued  wheel  surface,  the  maximum  wheel- soil  friction 
is  assumed  to  be  equal  to  the  internal  soil  friction. 


36 


The  value  of  the  constant 


«v»  *v»"!  W*  vrmVt  v*^*****,. 


t/  ,x  =  ^  •"■  ■n^^S  tan  $ 
( S=0)  s  +  m  ~ 


(15) 


■;  fits  best  in  both  cases  is  m  =  0.05  . 


Thus, 

s  =  3$  ;  t(&=0)  =  0.394  tan  $ 
s  =  14$  ;  ^5=0)  ~  0*774  tan  $ 

s  =  33$  ;  t(s=0)  =  °*^]L3  tan  ^ 

These  numerical  values,  although  determined  for  a  particular  point,  were 
considered  representative  of  the  total  zone  of  the  lunule 

*(8=0)  ”  H  (l6) 

Bulldozing  zone 

59*  In  the  bulldozing  zone,  all  soil  movements  occur  in  the  planes 
of  wheel  rotation;  there  is  no  lateral',  flow.  As  pointed  out  earlier,  the 
soil  in  the  immediate  neighborhood  of  the  wheel  within  the  bulldozing 
zone  is  in  a  perfect  plastic  condition.  The  rolling  surface  of  the  wheel 
is  assumed  to  be  a  rupture  surface,  and  the  t  value,  therefore,  is 

t^  =  tan  t  (17) 


It  is  known  that  along  the  center  line  of  the  contact  area,  the  bulldozing 
zone  extends  from  the  leading  edge  back  1 o  the  edge  of  the  lunule.^  It 
is  obvious,  however,  that  the  bulldozing  cannot  cover  the  entire  width  of 
the  wheel  everywhere.  Based  on  the  assumptions  already  made,  this  would 
lead  to  extremely  high  overall  tm  values,  for  both  low  and  high  slips.* 


*  If  it  were  assumed  that  the  bulldozing  zone  covered  the  whole  width  of 
the  wheel,  it  would  largely  predominate  over  the  lunule  zone  at  low  slip 
values.  Since  t  =  tan  is  assumed  for  the  bulldozing  zone,  the  tm 
value  of  the  total  wheel- soil  interface  would  be  close  to  tan  $  .  At 
high  slips,  the  lunule  zone  dominates.  Since  at  high  slips  the  t  value 
of  the  lunule  zone  is  near  the  maximum  tan  /  ,  the  overall  tm  value 
is  close  to  the  maximum. 
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60.  Therefore,  an  assumption  concerning  the  lateral  limits  of  the 
bulldozing  zone  was  necessary.  It  is  felt  that  the  increase  of  the  bull¬ 
dozing  effect  with  the  contact  angle  is  faster  than  a  linear  rate.  Hence, 
the  width  of  the  bulldozing  zone,  w_  ,  is  assumed  to  increase  with  the 
positional  angle  by  a  parabolic,  rather  than  a  linear,  function.  The 

simplest  expression  of  this  kind,  is  obtained  by  setting  w_  =  0  for 

B 

5=0,  and  w  =  b  for  5=1  (radian)  .  A  funnel-shaped  bulldozing 
zone  results  (fig.  19b),  the  equation  for  which  is 


w„  =  b6  (5  in  radians) 
B 


(18) 


Equation  18  was  chosen  for  its  mathematical  simplicity,  but  it  also  has 
physical  meaning,  since  5=1  radian  (=»57  deg)  is  roughly  the  maximum 
contact  angle  that  can  be  obtained  under  reasonable  operating  conditions. 

Ec  nation  18  states  that,  even  at  the  leading  edge  of  the  wheel,  the  bull¬ 
dozing  zone  does  not  cover  the  full  width,  b  ,  of  the  wheel  when  3  <  1 
radian  .  This  conclusion  is  justified  by  experimental  evidence,  because  at 
the  leading  edge  of  the  wheel,  the  soil  is  elevated  above  the  original 
level  (bow  wave)  end  tends  to  flow  laterally.  By  definition,  lateral  flow 
is  excluded  within  the  bulldozing  zone. 

61.  The  reason  that  the  origin  of  the  parabolas  was  put  at  the  bot¬ 
tom  dead  center  (5  =  0),  was  because  the  bulldozing  zone  at  its  real 

origin  (6  =  sa)  has  to  possess  a  certain  real  width;  otherwise,  it  could 

14 

not  have  bee1-  observed  by  Wong  and  Reece.  This  created  a  fictive  por¬ 
tion  of  the  *idozing  zone  within  the  lunuie  zone  (dashed  lines,  fig. 

19b),  which  was  not  considered  in  the  analysis. 

Lateral  flow  zone 

62.  The  two  remaining  portions  of  the  contact  area  (fig.  19b)  are 
termed  zones  of  lateral  flow  because  the  soil  particles  beneath  these 
zones  are  assumed  to  move  strictly  parallel  to  the  vertical  plane  through 
the  wheel  axle.  If  there  is  lateral  flow,  the  soil  obviously  is  in  a  per¬ 
fect  plastic  condition,  but  this  does  not  necessarily  imply  that  the  shear 
stress  in  any  direction  is  equal  to  its  maximum  value.  In  fact,  even  in 

a  perfect  plastic  state,  every  value  of  shear  stress  between  t  =  0  and 


36 


t  =  x  =o  ban  $  can  be  obtained,  depending  on  the  orientation  of  the 
max 

plane  considered.  The  shear  stress  is  zero  in  the  planes  of  the  principal 
stresses,  and  maximum  on  all  rupture  planes.  Since  there  were  no  reasons 
to  assume  that  the  rolling  surface  of  the  wheel  within  the  zones  of  lat¬ 
eral  flow  was  a  rupture  surface,  the  x/cr  ratio  may  yield  any  value  be¬ 
tween  zero  and  tan  p. 

63.  This  is  a  major  difference  from  the  bulldozing  zone,  where  the 
rolling  surface  was  intuitively  taken  as  a  rupture  surface.  Since  the 
pattern  of  rupture  surfaces  beneath  the  zones  of  lateral  flow  is  unknown, 
the  shear  stress  at  the  rolling  surface  cannot  be  assessed.  It  was  as¬ 
sumed  that  the  ratio  of  the  real  shear  stress  to  the  maximum  shear  stress 
at  the  rolling  surface  yields  a  constant  value  X  throughout  the  zone 

of  lateral  flow: 

r  =  Xx  =  Xcj  tan  $  (19) 

max  '  ‘ 

where  X  is  an  experimental  constant. 

64.  Another  fundamental  difference  between  the  zone  of  lateral  flow 
and  the  bulldozing  zone  is  that  the  real  shear  stress  vector  does  not  fall 
into  the  vertical  plane  of  travel,  but  has  a  significant  lateral  component 
due  to  soil  flow. 

65.  Since  only  the  tangential  component  in  the  vertical  plane  of 
travel  of  the  shear  stress  vector  contributes  to  the  pull  and  torque  char¬ 
acteristics  of  the  wheel,  the  effective  shear  stress  (i.e.  tangential  to 
the  wheel  perimeter)  is 

x  —  t  cos  7  =  X t  cos  7  =  Xc  tan  $  cos  7  (20) 

F  max  r  ‘ 

where  7  is  the  angle  of  inclination  of  the  t  vector  with  respect  to  the 
vertical  plane  of  travel. 

66.  In  other  words,  the  effective  shear  stress  within  the 

F 

zones  of  lateral  flow  is  reduced  with  respect  to  x  for  two  reasons: 

max 

(a)  the  x  vector  has  a  radial  component,  accounted  for  oy  the  coef- 

ficient  X  ;  (b)  the  x  vector  also  has  a  lateral  component,  accounted 

max 


for  by  the  cos  y  term.  At  present,  there  is  no  way  to  assess  the  nu¬ 
merical  value  of  \  ,  but  the  angle  y  might  be  roughly  evaluated  by 
the  following  reasoning. 

67.  A  local  velocity  vector  v^  can  be  attributed  to  every  point 
of  the  rolling  surface  (fig.  20a).  This  velocity  vector  is  perpendicular 
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Fig.  20.  Schematic  diagram  of  soil-wheel  velocity  vectors 


to  the  line  connecting  the  instantaneous  center  of  rotation  I  with  the 
point  P  in  question,  and  its  magnitude  is  as  times  the  length  of  the 
line  IP  . 

68.  The  velocity  vector  v  can  be  separated  into  its  tangential 

(v,  )  and  radial  (v  )  components,  which  are  (fig.  20) 
t  r 

v^  =  Rd)[l  -  (l  -  s)  cos  6]  (21) 


v^  =  Ros(l  -  s)  sin  & 


(22) 
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Only  the  wheel  displacements  due  to  vr  push  the  soil  away  from  under  the 

wheel,  while  tangential  wheel  displacement  is  a  pure  sliding  process  with 

respect  to  the  adjacent  soil.  The  direction  of  the  soil  movements  has 

already  been  assumed  to  be  purely  lateral.  It  is  further  hypothesized  that 

the  mean  horizontal  velocity  of  the  moving  soil  particles  v  is  equal  to 

s 

the  radial  wheel  velocity  that  causes  the  soil  to  move: 


v  =  v  =  Ru)(l  -  s)  sin  6 
s  r 


(23) 


Thus,  the  actual  soil-to-wheel  sliding  velocity  vgw  is  the  sum  of  the 
vectors  in  the  plane  of  the  rolling  surface  of  the  two  velocities  that 


contribute  to  it,  i.e.  v  and 
’  s 


vt  (fig.  20b).* 


V  =  V  +  V, 

sw  s  t 


Since  the  orientation  of  the  soil-to-wheel  shear  stress  vector  T  is  the 
same  as  that  of  the  soil-to-wheel  sliding  velocity  vector,  its  inclination 
with  respect  to  the  vertical  plane  of  travel  is  (fig.  20b) 


tan  y  = 


_  Vr  (l  -  s)  sin  5 
~  v,  ”  1  -  (1  -  s)  cos  5 


(24) 


Hence, 


cos  7  = 


1  -  (l  -  s)  cos  5 


Vl  +  tan2  7  ’ill  -  (1  -  s)  CCS  5]2  +  [(l  -  s)  sin  6]2 


(25) 


According  to  equations  20  and  25, 

T„ 


tF  = 


1  [1  -  (l  -  s)  cos  5] 


•jl  1  -  (l  -  s)  cos  6]2  +  [(l  -  s)  sin  5]2 


tan  $ 


(26) 


*  The  vector  v-^  is  drawn  in  the  opposite  sense  than  in  fig.  20a.  It 
represents  here  the  relative  velocity  of  the  soil  with  respect  to  the 


wheel  surface. 
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where  t_  is  the  effective  t  value  in  the  zone  of  lateral  flow. 

F 

69.  The  combined  width  of  both  zones  of  lateral  flow  (fig.  19b, 
equation  18)  can  be  expressed  by 

wp  =  b  -  vB  =  b(l  -  5)2  (27) 

The  length  of  the  two  lateral  zones  as  well  as  that  of  the  bulldozing  zone 
is  (fig.  19) 


ip  =  iB  =  (1  -  s)Ra 


(28) 


Plastic  and  elastic  regions 

70.  For  the  remainder  of  the  analysis,  the  bulldozing  zone  and  both 
lateral  zones  will  be  treated  jointly.  Together,  they  constitute  the 
principal  forward  or  plastic  region  because  the  adjacent  soil  is  in  a 
plastic  state.  Correspondingly,  the  principal  rearward  region,  which  is 
equivalent  to  the  lunule  zone,  is  the  elastic  region  (fig.  19). 

71.  From  equations  17,  18,  2 6,  and  27,  the  average  t  value  in  the 
plastic  region,  taken  over  the  wheel  width,  is 


t 

a 


Vb5 


t  = 
a 


_  (1  -  f)a  -  (1  -  s)  cos  5] 


+  6 


fa.  -  (1  >  s)  cos  5]2  +  [(l  -  s)  sin  5]2 


tan  0  (29) 


A  similar  equation  that  describes  the  variation  of  t  in  the  rearward 
elastic  zone  is  not  available.  It  is  assumed  that  equation  1 6,  which  is 
the  measurement  at  the  bottom  dead  center  (S  =  0),  can  be  considered  as 
the  mean  t/cj  value  over  the  total  elastic  zone. 

72.  Equations  16  and  29  are  represented  in  figs.  21  and  22  using  the 
following  numerical  values: 

tan  $  =  0.5  (reference  value) 

m  =  0.05  (experimental  value,  figs.  6  and  7) 
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.  22.  Theoretical  t 
tribution  at  soil- 
el  interface.  X  =  0.6 


X  =  0.4,  0.6 
s  =  3f>,  M>  33# 

where  X  is  an  experimental  coefficient;  the  numerical  values  of  0.4  and 
0.6  were  chosen  as  a  first  trial. 


Theoretical  and  Experimental  Results  Compared 
Using  New  Concept 

73-  The  theoretical  t  values  derived  using  the  proposed  approach 
were  compared  with  the  experimental  M/RW  values  using  equations  11a  and 
12b  in  the  same  manner  as  was  done  previously  usi~g  Sela's  theory.  Ex¬ 
amples  of  the  graphs  from  which  the  overall  mean  t  values  were  deter¬ 
mined  are  given  in  fig.  23.  These  plots  also  show  that  for  higher  con¬ 
tact  angles,  the  variation  of  the  t  ratio  generally  corresponds  with 
Sela's  experimental  data,  i.e.,  high  t  values  at  both  ends  of  the  con¬ 
tact  length,  minimum  value  near  the  middle.  For  small  contact  angles 
(a  +  p  <  0.5),  which  were  not  investigated  by  Sela,  the  t-variation  tended 
to  increase  from  the  leading  edge  toward  the  rear  end  of  the  contact 
length. 

74.  The  variations  of  t  with  the  contact  angle  (a  +  p)  are  il¬ 
lustrated  in  fig.  24.*  The  theoretical  t  curves  drawn  with  X  =  0.^ 
give  somewhat  lower  values  and  exhibit  a  slightly  steeper  increase  in 
their  ascending  portion  than  the  curves  of  X  =  0.6  .  When  these  theoreti¬ 
cal  curves  are  compared  with  the  experimental  results  (figs.  Il-l4),  it 
is  apparent  that  the  two  sets  of  curves  vary  in  almost  the  same  manner, 
as  required  by  equation  12b.  In  addition  to  the  complete  agreement  of  the 
five  qualitative  features,  a  rather  subtle  feature  of  the  experimental 
M/jRW  curves  is  reflected  by  the  t  curves:  for  s  =  14  percent  ,  the 
minimum  M/RW  value  occurs  at  a  higher  contact  angle  than  it  does  for  the 
two  other  curves,  except  in  fig.  11. 


*  To  avoid  confusion,  it  should  be  noted  t ne*  fig.  23  represents  the  varia¬ 
tion  of  t  as  a  function  of  the  position  angle  6  for  a  given  contact 
angle  (a  +  p);  fig,  24  gives  the  variation  of  t  as  a  function  of  the 
varying  contact  angle  (a  +  p). 
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Fig.  23.  Examples  of  graphs  used  to  determine  overall 
mean  t  values 


75*  As  a  second  check,  equation  11a  was  applied  to  the  proposed 
theory.  It  was  found  that,  for  tests  run  with  a  cone  index  of  60,  the 
X  =  0.4  curves  fit  better;  for  a  cone  index  of  30,  the  X  =  0.6  curves 
have  a  better  fit.  Using  the  appropriate  X  values  and  correcting  the 
value  of  tan  $  as  previously  explained,  the  experimental  M/RW  curves 
are  plotted  in  figs.  25  and  2 6,  together  with  the  theoretical  t  curves 
and  the  +25  percent  error  limits.  With  the  exception  of  some  portions 
of  the  3  percent  slip  curves,  the  M/RW  curves  fall  within  the  theoreti¬ 
cal  bounds.  Only  a  small  percentage  of  the  total  length  of  all  curves 
lies  beyond  the  limits.*  This  small  error  percentage  can  be  reduced  even 
further  by  a  more  appropriate  choice  of  the  X  value  (see  paragraph  79)- 
There  is  good  agreement  between  the  theoretical  values  and  the  experi¬ 
mental  findings.  Thus,  the  basic  idea  leads  to  a  reasonably  valid  approxi¬ 
mation,  at  least  for  the  range  of  conditions  tested. 

Comparison  of  Hew  Concept  and  Sela's  Approach 

76.  Many  of  the  basic  assumptions  of  the  proposed  approach  are  01- 

17 

posite  to  those  made  by  Sela.  Sela's  theory  is  essentially  two- 
dimensional  (lateral  flow  is  neglected);  the  new  concept  is  based,  to  a 
great  extent,  on  consideration  of  the  lateral  flow.  Furthermore,  an  at¬ 
tempt  was  made  to  take  into  account  the  effect  of  the  lunule,  the  existence 
of  which  was  unknown  at  the  time  Sela's  theory  was  elaborated. 

77.  The  most  Significant  difference,  however,  is  the  rejection  of 
the  idea  that  the  shear  stress  depends,  among  other  factors,  on  the  "dis¬ 
tance  from  the  point  in  question  to  the  leading  edge  of  the  contact 
area."  When  this  is  assumed,  the  maximum  x/cr  ratio  is  never  reached, 
even  at  100  percent  slip,  which  contradicts  common  sense.  Under  the  new 
concept,  equations  l4  and  15  indicate  that  for  100  percent  slip,  x/cr 

=  tan  $  over  the  total  area  of  the  wheel-soil  interface.  However,  the 
theoretical  concept  developed  herein  is  not  proposed  as  a  fi..al  answer 

*  As  explained  for  Sela's  results,  this  error  percentage  would  not  change 
significantly  if  the  experimental  M/RW  curves  were  shifted  slightly  to 
the  right  to  compensate  for  the  faulty  measurement  of  contact  length. 
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G  =  2.72  N/ CM 


Fig.  25.  Corrected  tjn  curves  with  +2%  error  limits 
superimposed  on  corresponding  m/RW  curves,  X  =  0.6 


to  the  problem  of  the  soil-wheel  interaction.  Some  very  crude  assumptions 
had  to  be  made  that  will  require  revision  as  additional  data  are  accumu¬ 
lated.  Rather,  this  concept  is  intended  to  demonstrate  that  an  entirely 
different  approech  is  possible. 

78.  The  most  controversial  assumption  of  the  new  concept  is  that 
of  a  constant  X  value  valid  over  the  whole  zone  of  lateral  flow.  It 
should  he  remembered  her  •  that  the  parameter  X  was  introduced  in  an  at¬ 
tempt  to  comprehend  the  influence  of  a  complex  physical  phenomenon — the 
rupture  pattern  beneath  the  zone  of  lateral  flow — by  a  single  constant. 
This  hypothesis,  hovever,  makes  the  approach  very  flexible  because  the 
parameter  X  is  an  experimental  coefficient  and  can  he  chosen  so  that 

the  theoretical  t  curves  approach  the  M/RW  curves  even  closer  than 
m 

in  figs.  25  and  2 6.  The  form  of  the  t  curves  is  not  changed  signifi- 

m 

cantly  by  a  varying  X  (fig.  24),  hut  the  position  is.  As  X  increases, 
the  t  curve  shifts  upward.  Thus,  the  somewhat  unexpected  lower  M/RW 
values  at  higher  cone  index  ratings  (compare  figs.  13  and  11;  l4  and  12) 
could  effectively  be  dealt  with  by  associating  the  lower  of  the  two  X 
values  with  the  higher  soil  strength  and  vice  versa. 

79-  Even  if  further  investigations  verify  this  concept  of  a  con¬ 
stant  X  value  over  the  whole  zone  of  lateral  flow,  it  is  probable  that 
X  varies  not  only  with  soil  strength,  but  also  with  slip,  wheel  geometry, 
and  load.  The  determination  of  X  can  be  refined  using  available  test 
data.  For  example,  the  narrow  wheel  appears  to  develop  a  higher  M/RW 
ratio  than  the  wide  wheel  (compare  figs.  11  and  12  with  figs.  13  and  l4). 
Hence,  X  decreases  as  wheel  width  increases.  By  similar  arrangement,  an 
almost  perfect  quantitative  coincidence  of  t  and  M/RW  curves  could  be 
obtained  for  the  3  percent  slip  condition,  instead  of  the  rather  poor 
agreement  between  curves  for  the  chosen  values  of  X  shown  in  figs.  25 
and  26. 

80.  Such  considerations,  however,  are  somewhat  premature.  This 
concept  should  be  thoroughly  investigated,  and  hopefully,  the  basic  as¬ 
sumptions  refined,  before  a  definite  pattern  of  numerical  values  is 
established. 


50 


PART  VI:  CONCLUSIONS  AND  RB3CM4KNDATIOBS 


Conclusions 

8l.  The  following  conclusions  are  believed,  warranted  by  the  findings 
of  this  report: 

a.  The  shear  stress-displacement  relation  for  soils  (equa¬ 
tion  2  or  similar  ones)  is  inadequate  for  describing  the 
soil-wheel  interaction,  or  for  predicting  wheel  performance. 

b.  The  mean  shear-to-normal  stress  ratio  tffl  at  the  wheel- 
soil  interface  is  closely  related  to  the  M/few  ratio 
(equations  11a  and  12b),  and  can  be  used  to  experimentally 
determine  the  validity  of  any  t/o  theory  more  accurately 
than  can  be  done  using  direct  measurements  of  the  tangen¬ 
tial  and  normal  stresses. 

c.  A  new  concept,  based  on  the  -real  soil  rupture  pattern  be¬ 
neath  rigid  wheel.' ,  provides  a  sound  approach  for  deter¬ 
mining  the  t/c  ratio. 


Recommendations 


82.  It  is  recommended  that: 

a.  Experimental  and  theoretical  research  be  conducted  to  study 
the  flow  pattern  beneath  rigid  wheels,  especially  the 
three-dimensional  features  of  the  pattern.  By  using  a  rigid 
wheel  with  a  transparent  rolling  surface,  the  flow  pattern 
imaediately  beneath  the  wheel  could  he  photographed  from 
inside  the  wheel.  Thus,  the  assumptions  concerning  the 
different  zones  at  the  wheel-soil  interface  could  he  checked 
directly  and  refined. 

h.  Investigations  of  the  t/o  ratio,  similar  to  those  reported 
here,  he  made  using  cohesive  soil. 

c.  In  conjunction  with  the  research  on  the  t/o  ratio,  sys¬ 
tematic  measurements  he  made  of  the  variations  of  the  nor¬ 
mal  stresses  within  the  wheel-soil  interface,  and  these 
variations  be  analytically  expressed.  Once  the  normal 
stress  and  the  t/o  ratio  can  be  predicted,  the  theoreti¬ 
cs!  problems  of  soil-wheel  mechanics  for  rigid  wheels  can 
be  considered  solved.  This  solution  appears  to  he  pre¬ 
requisite  to  considering  the  more  complex  problems  of  de¬ 
formable  tires  on  soft  soils  from  other  than  a  purely  em¬ 
pirical  point  of  view. 
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\  i»  amtract 

Current  theoretical  concepts  of  soil-running  gear  interaction  are  based  on  empirical 
pressure-sinkage  and  shear  stress-displacement  relations.  However,  tests  from  which 
these  relations  were  obtained  present  a  poor  analogy  to  the  soil-running  gear  interac¬ 
tion.  In  particular,  use  of  the  shear  stress-displacement  relation  for  analytical 
evaluation  of  the  shear  stresses  at  the  contact  surface  of  rigid  wheels  on  sand  was 
found  misleading,  primarily  recause  equations  obtained  by  the  bevameter  or  dragged 
plate  tests  describe  soil  benavior  in  an  intermediate  state  between  the  elastic  and 
plastic  states.  Recent  publications  show  that  part  of  the  soil  close  to  a  powered 
wheel  is  in  an  actual  plastic  flow  state  and  part  is  in  a  quasielastic  state.  Thus, 
the  soil-bevameter  interaction  is  fundamentally  different  than  the  soil-wheel  interac¬ 
tion.  Sela's  theory  of  the  relation  between  a  rigid  wheel  and  dry  sand  is  based  ex¬ 
clusively  on  the  shear  stress-displacement  concept,  and  thus  provided  an  excellent 
means  of  checking  the  concept.  Using  a  simple  approximate  relation  between  the  M/RW 
ratio  and  the  mean  shear-to-ncrrmal  stress  ratio  (t/o)  taken  over  the  total  contact 
surface,  the  theory  was  checked  by  means  of  constant-slip  rigid  wheel  tests.  Results 
strongly  supported  the  conclusion  that  shear  stress-displacement  relations  are  ir¬ 
relevant  to  soil-wheel  mechanics .  A  new  theory  was  developed  that  attempts  to  assess 
the  variation  of  t/o  along  the  soil-wheel  interface,  without  referring  to  stress- 
displacement  relations. .-It  uses  the  soil  rupture  pattern  beneath  rigid  wheels  on  sand 
to  subdivide  the  total  interface  into  three  zones  of  different  soil  behavior.  This 
theory  agrees  better  with  teat  data  than  other  theories. 
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Errata  Sheet 
Wo.  1 

CONTRIBUTION  TO  THE  MECHANICS  OF 
RIOS)  WHEELS  Of  SAND 

Technical  Report  M-68-2 

May  1968 

1.  On  page  7,  line  3,  change  "critical"  to  "rupture." 

2.  On  page  8,  in  title  of  fig.  2,  raised  number  should  be  "13,"  not  "7." 

3.  On  page  18,  in  fig.  9  tinder  "t-DIS  RIBUTICNS,"  change  "z"  to  "n"  in 
all  cases. 

4.  On  page  20,  line  2  of  paragraph  40,  insert  "the"  between  "from"  and  "t  ." 

5.  On  page  21,  line  7  below  equation  12,  change  '(fig.  9)"  to  "(fig.  10)." 

?>.  On  page  44,  last  line,  change  "fig.  11"  to  "fig.  14." 
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